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ABSTRACT 

(N 

>- ■ We present initial results from the first optical spectroscopic survey of high latitude 

^ . Galactic cirrus clouds. The observed shape of the cirrus spectrum does not agree with 

Q i that of scattered ambient Galactic starlight. This mismatch can be explained by the 

CN I presence of Extended Red Emission (ERE) in the diffuse interstellar medium, as found 

in many other astronomical objects, probably caused by photoluminescence of hydro- 
Q\ . carbons. The integrated ERE intensity, /ere ~ 1-2 x 10~^ergs^^ cm~^sr~^, is roughly a 

third of the scattered light intensity, consistent with recent color measurements of diffuse 
Ph| Galactic light. The peak of the cirrus ERE (Aq ~ 6000 A) is shifted towards short (bluer) 

Q ■ wavelengths compared to the ERE in sources excited by intense ultraviolet radiation, 

such as Hll regions (Aq ~ 8000 A); such a trend is seen in laboratory experiments on 



i-i 



j^ i hydrogenated amorphous carbon films. 

>: 

k> , Subject headings: Galaxy: solar neighborhood-ISM: clouds-ISM: dust, extinction-ISM: 

H ' molecules-ISM: reflection nebulae 
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1. Introduction 

A diffuse component to the Galactic interstellar 
medium (ISM) was noticed over forty years ago (de 
Vaucouleurs 1955, 1960) in the form of scattered 
starlight from dust clouds (Sandage 1976). It was 
the IRAS all-sky survey though that first revealed cir- 
rus as a ubiquitous component of the far-infrared sky 
(Low et al. 1984), demonstrating that these tenuous 
clouds of warm dust are representative of the diffuse 
ISM in our Galaxy (Beichman 1987; Boulanger 1994) 
and probably in all late-type spirals. Low et al. (1984) 
noted the positional correlation of the cirrus emission 
with Hi clouds mapped by Heiles (1975); later optical 
studies of cirrus (de Vries & le Poole 1985; Paley 1990; 
Stark 1993) revealed a good correspondence between 
optical and infrared cirrus. 

Since the IRAS mission, analyses of the cirrus have 
shown that the far-infrared emission from high lati- 
tude cirrus arises from an ensemble of silicate and 
graphite dust grains and complex molecules heated 
by the ambient interstellar radiation field (ISRF) 
(Low et al. 1984; Mezger et al. 1982). The color 
temperature of the cirrus emission at 12 and 25/zm, 
where visible, indicates thermal fluctuations of very 
small grains and/or polycyclic aromatic hydrocar- 
bons (PAHs) within clouds (cf. Puget et al. 1985; 
Boulanger et al. 1985, 1988; Heiles et al. 1988; Guha- 
thakurta & Draine 1989). There is at least an order of 
magnitude variation in the infrared colors from cloud 
to cloud, indicating considerable compositional vari- 
ation. 

In addition to scattered and thermally reprocessed 
radiation, Extended Red Emission (ERE), a broad 
emission feature of ~ 1000 A width centered between 
6000 A and 8000 A, has been detected in the spectra 
of a wide variety of Galactic objects: reflection nebu- 
lae (Schmidt et al. 1980; Witt & Schild 1988), a high 
latitude Galactic dark cloud (Mattila 1979; Chlewicki 
& Laureijs 1987), planetary nebulae (PN) (Furton & 
Witt 1990, 1992), Hii regions (Perrin & Sivan 1992; 
Sivan & Perrin 1993), and a nova (Scott et al. 1994). 
Recently, ERE has also been detected outside the 
Galaxy, in the halo of M82 (Perrin et al. 1995). 

Possible carriers of the ERE include hydrogenated 
amorphous carbon (HAG — Duley 1985; Duley & 
Wilhams 1988, 1990), PAHs (d'Hendecourt et al. 
1986), and Geo (Webster 1993). There is no obvious 
"best" candidate. PAHs have been identified with 
the so-called unidentified infrared emission (UIR) fea- 



tures and indeed, in the reflection nebula NGC 2023, 
emission at 3.3 /im is correlated spatially with ERE 
(Burton et al. 1988; Witt & Malin 1989). However, 
Perrin & Sivan (1992) do not find such a correlation in 
the Orion nebula; moreover, 3.3 /im emission has been 
detected in PN with no ERE (Furton & Witt 1992). 
Laboratory spectra of Ceo compare well with the ob- 
served ERE in the PN NGC 2327 (Webster 1993). 
The peak and shape of the ERE feature though vary 
for different objects, and it is unclear whether Ceo can 
reproduce this variability. Whether ERE is present in 
all dust or only in discrete dusty objects remains, to 
some degree, an open question. The study of ERE in 
the warm diffuse ISM has important implications for 
the chemistry of and energy balance in this phase. 

There is indirect evidence of ERE in the Galaxy's 
warm diffuse interstellar medium. Unusually red 
broadband colors {B — R, R — I) have been measured 
for cirrus clouds (Guhathakurta & Tyson 1989; Paley 
et al. 1991; Guhathakurta & Gutri 1994), suggesting 
the presence of an additional component (over scat- 
tered light) in the R- and /-bands. Recent work by 
Gordon et al. (1997) reached a similar conclusion us- 
ing a very different technique: they decomposed Pio- 
neer 10 & 11 blue and red surface brightness maps into 
contributions of Galactic starlight and diffuse Galac- 
tic light. The advantage of this approach is the wide 
field of view, while the obvious disadvantage is poor 
angular resolution. All of these studies lack the spec- 
tral resolution needed to make an unambiguous detec- 
tion of the ERE bump. Their interpretation is further 
complicated by the lack of detailed knowledge of the 
color of incident starlight and the amount of internal 
reddening within clouds. 

This paper presents the first direct spectroscopy 
of diffuse, high latitude cirrus clouds illuminated by 
ambient Galactic starlight. The targets of study have 
low optical surface brightness, only ~ 1% of the night 
sky brightness (they are not bright enough to be clas- 
sified as reflection nebulae), but the experiment is 
made possible by the availability of high efficiency 
modern spectrographs. These cirrus clouds have op- 
tical depths substantially below unity (typically a few 
tenths) making them qualitatively different from the 
previously studied dark cloud, L1780 (Mattila 1979), 
with Ab ^ 2-5 mag. Our spectroscopic survey covers 
the relatively tenuous outer edges of cirrus clouds, but 
even their densest central portions show none of the 
usual signs of high opacity: reduced surface density 
of background stars or increased IRAS lOO/^m/60/im 



ratio. Sec. 2 contains a description of the observa- 
tions and data reduction techniques. In Sec. 3, the 
spectra are analyzed in the context of scattering of 
starhght and photoluminescence from hydrocarbons. 
The main findings of this paper are summarized in 
Sec. 4. 

2. Observations &: Data Reduction 

Our targets of study arc compact, relatively iso- 
lated, high latitude cirrus clouds selected from IRAS 
100 /Ltm maps or drawn from the catalog of Lynds 
Bright Nebulae (Lynds 1965) which is, in turn, de- 
rived from POSS plates. Clouds with sharp edges 
are ideal for making (differential) measurements of 
clouds relative to the adjacent sky background. The 
Lick Observatory Shane 3-meter telescope and KAST 
spectrograph were used to observe three cirrus clouds 
- Witch (10'^44'"42% -h83°43 f7), LBN8 (15'^50'"3% 
-5°50.'2), and Draco (16'^47'"34% -F6G°17 .'0) ~ dur- 
ing 2 nights in April 1997. The KAST spectrograph's 
300 fines mm^^ grating and Reticon 1200x400 CCD 
(of wfiicfi 1200x200 are read out) yield a wavelength 
range of AA = 4090A-9590A at a dispersion of 
4.6 A pixel^^ with 27 /im pixels and a scale of about 
'.' 8 pixer 1 . The slit width and length of 3" x 145" re- 
sult in an overall resolution of 18. 4A (FWHM). The 
observations consist of 20 min exposures x 15 point- 
ings for the Witch and 7 pointings for each of LBN 8 
and Draco, with each pointing chosen to have the slit 
straddle a sfiarp edge in the cirrus. Dome flat spectra 
were obtained after every change of pointing and/or 
of slit position angle. Tfie spectropfiotomctric stan- 
dard stars Feige 34 and Fcige 98 were observed several 
times during the observing run. Observing conditions 
were initially photometric but degraded during the 
course of the first night; the second night was not 
photometric. The seeing FWHM was ~ 2" during 
the observations. 

As both zero level and dark current are negligible, 
we do not apply any overscan, bias, or dark current 
subtractions. For each spectrum, a flat-field image is 
constructed by averaging the best matched dome flats 
(nearest in terms of telescope pointing and slit posi- 
tion angle). Due to flexure of the instrument however, 
the fringe patterns in spectra and flat fields do not 
correspond exactly. To correct for this effect, we al- 
low for a small shift of the flat fields along the disper- 
sion axis. The optimal shift, typically ~ 1-2 pixels, is 
determined by minimizing the rms brightness varia- 



tion along the slit of several bright night sky emission 
lines in the fiat-fielded images. The spectra are then 
corrected for geometric distortion in the spectrograph 
optics by rectifying the night sky lines. 
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Fig. 1. — Spectra of three cirrus clouds, in order of 
decreasing surface brightness. Representative error- 
bars are shown, spanning the smallest and largest er- 
rors at the blue and red ends of the spectrum. The 
solid lines are scattered synthetic E-type galaxy spec- 
tra (Poggianti 1997) normafized to the cirrus spectra 
at ~ 4500 A. Normalizing the scattered light model 
at 5000 A (dot-dashed line in first panel) reduces the 
inferred amount of ERE but results in a deficit for 
A < 4800 A; the dotted line in the second panel is 
based on an Sc-typc ambient spectrum; the long- 
dashed line in the third panel corresponds to twice 
the nominal optical depth. The lowest panel shows 
the night sky spectrum. 

The low surface brightness of the cirrus {B ~ 
26magarcsec~^, R ~ 24magarcsec^^) makes this a 
challenging project, and requires very accurate sub- 
traction of the night sky spectrum, especially its 
bright emission lines. For a few of the spectra the 
task is facilitated by the fact that the long slit ex- 
tends to cirrus-free parts of the sky on both sides 
of the cloud. The other cases are more complicated, 
and we use a multi-step iterative process: [1] A rel- 



atively emission line free part of the night sky spec- 
trum (6350 A-7515 A) is collapsed in the dispersion 
direction (combining several columns) to determine 
the exact location of the cirrus edge (^cdgo ~ 0.5 ?max), 
and the brightness profile of the cirrus along the slit 
length, I{1). [2] The cirrus free section [0, ledge] is 
used to make a preliminary measurement of the sky 
spectrum, which is then subtracted from each row of 
the full 2D spectrum. [3] This extrapolation of the 
night sky spectrum beyond / > /edge has large un- 
certainties, however, primarily because of imperfect 
rectification. We isolate a small (A/ ~ 10") section of 
the (preliminary) sky-subtracted spectrum at the end 
of the slit opposite the cirrus-free end, [0.9 /max, /max], 
and make a low-order polynomial fit, /li", to the cir- 
rus spectrum (relying on the assumption that it is 
fairly smooth in A). [4] This smooth fit, I^"'^, is used 
to subtract off the contribution of the cirrus from the 
[0.9 /max, ^max] scction of thc "original" 2D spectrum 
(i.e., the pre-sky subtraction, flat-fielded and rectified 
spectrum). The final night sky spectrum is then in- 
terpolated from two extreme sections of this residual 
2D spectrum: [0, /edge] and [0.9 /max, /max]- This fi- 
nal night sky spectrum is subtracted from each row 
of the 2D spectrum. 

Next, each sky-subtracted spectrum is collapsed 
into a median cirrus spectrum over the appropri- 
ate section of the slit length: [/edge, 0.9 /max]- Me- 
dian spectra from individual exposures of a given 
cloud are then combined using a noise-weighted av- 
erage, with the noise determined over the "line-free" 
6350 A-7515 A range. The rms scatter in each A pixel 
amongst individual spectra is used to mask poorly 
subtracted night sky emission lines, especially at the 
red end where residual fringes are a problem. The 
resulting masked spectrum is then boxcar-smoothed 
in A, ignoring masked pixels, to a resolution of 97 A 
in order to enhance the signal-to-noise ratio. 

Wavelength calibration is done using night sky 
lines (Osterbrock & Martel 1992) in a star- and 
cirrus-free part of one of the spectra. Initial flux 
calibration is based on observations of the Feige 34 
spectrophotometric standard star. Second order over- 
lap, combined with the high flux of Feige 34 for 
A < 4000 A, leads to calibration errors for A > 7500 A. 
To correct this, we use the initially calibrated spectra 
of a number of field stars that the long slit happened 
to intersect. Most of these stars are easily identifiable 
as K and M dwarfs, as their observed spectra closely 
resemble empirical stellar library spectra (Fioc & 



Rocca-Volmerange 1997) for A < 7500 A. This enables 
us to determine a residual flux calibration correction 
for A > 7500 A (averaged over all the field stars) which 
is then applied to the cirrus spectra. 

3. Evidence for Extended Red Emission 

Figure 1 shows the calibrated spectra of all three 
cirrus clouds. We also plot a night sky spectrum to 
illustrate the strength of the night sky lines, partic- 
ularly for A > 7500 A. Since no particularly bright 
star is seen in the vicinity of the clouds, it is rea- 
sonable to assume that the cirrus is illuminated by 
ambient Galactic starlight. This assumption is jus- 
tified because: (1) the spectral shapes of clouds in 
different locations are roughly similar, and (2) our 
measurement represents an average over multiple slit 
positions within each cloud, effectively averaging out 
the illumination effects of individual bright stars. 

We compare the observed spectral shape of the 
cirrus to a simple scattering model. The spectrum 
of the incident radiation field is represented by syn- 
thetic galactic spectra of various Hubble types (Pog- 
gianti 1997). The optical depth estimate is based 
on a standard conversion from IRAS 100 /im bright- 
ness (Savage & Mathis 1979; Boulanger & Perault 
1988): Av = 0.053 (S'ioo/1 MJysr-^). This yields 
values of Ay = 0.27 for the Witch and Draco and 
Ay = 0.64 for LBN 8. In Fig. 1 we also plot the scat- 
tered starlight model normalized to the cirrus emis- 
sion at ^ 4500 A, based on an incident E-type galaxy 
spectrum. The grain albedo is taken to be constant 
in \ at uj — 0.55 (Draine & Lee 1984), while the ex- 
tinction law of Cardelli et al. (1989) is adopted. 

It is clear that the overall shape of the cirrus spec- 
trum is not fit by scattered starlight alone. In partic- 
ular, the cirrus spectra have too much tilt/curvature 
in the region 4000A-7000A to be fit by any plausi- 
ble mix of stars (spectra of galaxy types later than 
E produce even worse fits). Scaling the model spec- 
trum to match the red part of the cirrus spectrum 
(A ^ 5000 A) results in a deficit for A <, 4800 A 
(cf. top panel of Fig. 1); we have instead chosen to 
scale the model spectrum to the blue side of the cir- 
rus as the resulting broad region with excess emis- 
sion (ERE) has a natural explanation. It is not quite 
clear whether this excess emission extends beyond 
A ^ 7500 A where measurement/calibration uncer- 
tainties are large. Incidentally, the shape of the cir- 
rus spectrum agrees with previous broadband obser- 



vations of cirrus clouds (cf. Guhathakurta & Tyson 
1989), which showed both B — R and R — I colors to 
be too red to be explained by scattered Galactic light. 

Note, the 6*100 values adopted are conservative in 
that they have not had any large scale background 
subtracted (due to foreground zodiacal light, for ex- 
ample); for LBN 8 in particular, the actual brightness 
of the cloud, and the inferred Ay ^ could be as much 
as a factor of 2 lower. The conversion factor of 0.053 
used to derive Ay from 6*100 is well within the range 
found by Stark (1995) for Galactic dust clouds, based 
on photometry and counts of background stars. As an 
extreme example, we adopt a conversion factor of 0.1 
for the densest cloud in our sample, LBN 8. The re- 
sulting scattered light spectrum (long-dashed line in 
third panel of Fig. 1) is understandably somewhat (in- 
ternally) reddened relative to that derived using the 
nominal optical depth (solid line) , but the conclusions 
about ERE in the cirrus are essentially unchanged. 
The effect of doubling the nominal optical depth is 
negligible for the other two clouds. Multiple scatter- 
ing is ignored because Ty^ < 0.6 even for the densest 
cloud (LBN 8) using an extreme 6100 -^ Ay conver- 
sion ratio, whereas such events are only important if 

We also check whether the observed flux density 
of the cirrus is consistent with an independent esti- 
mate of the strength of the ISRF in the solar neigh- 
borhood (Mathis et al. (1983): 47r Ja(7000A) = 1.5 
and 3.1 ergs~^cm~^/im~^, at galactocentric distances 
of 10 and 8 kpc, respectively. Applying the above 
dust scattering parameters to the actual cirrus spec- 
tra we derive values for the incident ISRF of 1- 
2ergs^^cm^^/im^^. The agreement with Mathis et 
al.'s estimate is remarkable when one considers that: 
(1) it is difficult to measure the effective optical depth 
of the cirrus averaged over the exact area represented 
by the spectrum, given the poor angular resolution 
of the IRAS 100 /im maps; and (2) the "isotropic" 
ISRF strengths derived from the cirrus are likely to 
be biased low if we are primarily seeing back scat- 
tered radiation from preferentially forward scattering 
grains in clouds that are mostly illuminated on one 
side by the Galactic disk. 

The absence of a detectable scattered Ha emission 
line in the cirrus spectrum suggests that the radiation 
incident on the clouds has the spectrum of an early 
type galaxy, E or Sa. This is understandable if the cir- 
rus is local (heliocentric distance ^ 300 pc), in which 
case it is illuminated mainly by a part of the disk 



between spiral arms (which are separated by roughly 
2-3 kpc) . The inter-arm region is dominated by stars 
that are somewhat older than the global average over 
the Galaxy, the latter including a contribution from 
young stars in star forming/Hii regions. This expla- 
nation for the lack of Ha is indirectly corroborated by 
independent evidence that the cirrus ERE is excited 
by a UV-poor radiation field (see below). Moreover, 
the presence of ERE tends to dilute the reflected Ha 
emission line, if any, in the cirrus spectra. For com- 
pleteness, we also present the scattered light spectrum 
for an ambient Sc-type radiation field (dotted line in 
second panel of Fig. 1); it is a worse fit to the data 
than scattered light based on an E-type spectrum. 
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Fig. 2. — The difference between the cirrus spectrum 
and the best fit (at ~ 4500 A) scattered light model, 
averaged over the three clouds (bold line) and for 
LBN 8 and Witch (upper and lower thin solid lines). 
The dashed vertical line shows the approximate cen- 
tral wavelength of the ERE. The arrows show the 
peak wavelengths of the ERE in a number of astro- 
nomical objects; a reflection nebula (RN), a dark neb- 
ula (DN), a planetary nebula (PN) and Hii regions. 

In Fig. 2 we plot the average difference between 
the cirrus spectrum and the best fit (at ^ 4500 A) 



scattered light model, over the 4000A-7500A range 
within which flux measurements are reliable. The thin 
solid lines indicate the ERE for LBN 8 and Witch, the 
clouds with the strongest and weakest excess, respec- 
tively. These also provide some measure of the un- 
certainty in our estimate of the average cirrus ERE 
spectrum. 

Furton & Witt (1993) have conducted laboratory 
experiments on HAC and AC (pure amorphous car- 
bon) films to address astronomical observations of 
ERE. They find that UV radiation can either dimin- 
ish or enhance the efficiency of photoluminescence 
(PL) in hydrogenated carbon solids, depending on the 
availability of free atomic hydrogen and on the rela- 
tive H saturation of the solid. UV-annealing and sub- 
sequent hydrogenation of HAC dust causes a broad- 
ening and shift of the ERE peak towards the red. Fur- 
ton & Witt note that the observed shifting of the peak 
ERE emission to longer wavelengths in objects with 
intense UV radiation, such as Hii regions and plane- 
tary nebulae, must indeed be a consequence of this. 
Comparing the shape of the cirrus ERE (Fig. 2) with 
laboratory PL yield (Fig. 1 of Furton & Witt 1993), 
we see excellent agreement in both shape and width 
of the feature. We also indicate in Fig. 2 the approxi- 
mate central wavelengths Aq of the ERE in a number 
of astronomical objects: the Hii regions NGC 7635 
(Sivan & Perrin 1993) and the Orion Nebula (Perrin 
& Sivan 1992), the planetary nebula NGC 7027 (Fur- 
ton & Witt 1990), the dark nebula L1780 (Mattila 
1979; Chlewicki & Laureijs 1987), and the reflection 
nebula called the "Red Rectangle" around HD 44179 
(Schmidt et al. 1980). In agreement with Furton & 
Witt's laboratory experiments, the peak of the cirrus 
ERE, Ao ^ 6000 A, appears to be located at bluer 
wavelengths than the ERE emission seen in environ- 
ments with higher UV radiation. 

The integrated ERE intensity of the cirrus is 1.2 x 
10~^ ergs~^ cm~^ sr^^ and the ratio of ERE to scat- 
tered intensity is ^ 0.3, in good agreement with the 
findings of Gordon et al. (1997): /ere ~ lO^^ergs^^ 
cm"2 sr^i and Iere/Iscii == 0.05-2 in the diffuse ISM. 
For other objects Ieke/Isch ranges from 0.01-0.68 in 
reflection nebulae to 0.2-0.6 in Hii regions. Applying 
the same conversions as Gordon et al. (1997), we find 
the cirrus has an ERE photon efficiency of ~ 10%. 

The shape of inferred cirrus ERE spectrum looks 
remarkably similar to laboratory ERE spectra and to 
ERE in other astrophysical objects, and there is a 
perfectly reasonable explanation for the shift of the 



cirrus ERE peak to the blue relative to objects with 
hotter illuminating sources. Nevertheless, the cirrus 
ERE spectrum derived from this initial study should 
be treated as preliminary (pending more thorough 
follow-up observations) in light of the uncertainties 
in the data and in the scattered fight model discussed 
above. 

4. Summary 

We present long slit spectra straddling sharp edges 
in three high latitude Galactic cirrus clouds. The cir- 
rus is clearly detected in each case. The shape of the 
cirrus spectrum does not match the shape of scattered 
ambient interstellar radiation, modeled using syn- 
thetic galaxy spectra. Assuming the blue end of the 
cirrus spectrum is dominated by scattered light, there 
is an excess emission in the range 5000A-7000A, 
in good agreement with the Extended Red Emission 
(ERE) bump seen in other astronomical objects. The 
integrated ERE intensity {^ 10^'"' ergs~^ cm~^ sr~^) 
and ratio between ERE and scattered light (~ 0.3) 
agree well with values inferred from color measure- 
ments of the diffuse Galactic light (Gordon et al. 1997). 
The position of the peak ERE emission in cirrus 
(^ 6000 A), compared to objects such as Hii regions 
and planetary nebulae, is consistent with a relatively 
UV-poor ambient radiation field. 
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